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ABSTRACT 



We report on observations of the dust trail of comet 67P/Churyumov-Gera- 
simenko (CG) in visible light with the Wide Field Imager at the ESO/MPG 
2.2m telescope at 4.7 AU before aphelion, and at 24 /im with the MIPS instru- 
ment on board the Spitzer Space Telescope at 5.7 AU both before and after 
aphelion. The comet did not appear to be active during our observations. Our 
images probe large dust grains emitted from the comet that have a radiation 
pressure parameter /3 < 0.01. We compare our observations with simulated 
images generated with a dynamical model of the cometary dust and constrain 
the emission speeds, size distribution, production rate and geometric albedo 
of the dust. We achieve the best fit to our data with a differential size distri- 
bution exponent of -4.1, and emission speeds for a /3 = 0.01 particle of 25 m/s 
at perihelion and 2 m/s at 3 AU. The dust production rate in our model is 
on the order of 1000 kg/s at perihelion and 1 kg/s at 3 AU, and we require 
a dust geometric albedo between 0.022 and 0.044. The production rates of 
large (>10 /im) particles required to reproduce the brightness of the trail are 
sufficient to also account for the coma brightness observed while the comet 
was inside 3AU, and we infer that the cross-section in the coma of CG may 
be dominated by grains of the order of 60-600 /im. 



Keywords: COMETS, DUST; INFRARED OBSERVATIONS 
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1 Introduction 



A cometary trail is a narrow and elongated structure along the projected orbit 
of a comet. Trails are thought to consist of dust particles characterised by little 
sensitivity to solar radiation pressure and emitted at low speeds relative to 
the comet, such that they remain close to the orbit of the parent comet for 
many revolutions around the Sun. Such properties are commonly associated 
with mm- to cm-sized particles, and the trail is therefore interpreted as a 
continuous record of the emission history of such particles from the comet 



( ISykes and Walker! fjl992l )). The tail of a comet, by contrast, consists of smaller 
dust particles, which are more sensitive to radiation pressure and are therefore 
removed from the cometary environment on timescales of weeks to months. 
The tail roughly points away from the nucleus in the direction opposite to the 
Sun, while the trail is found along the orbit of the comet. Since the transition 
from "small" to "large" particles is likely a smooth one, the distinction between 
tail and trail is not always clear, in particular when the observer is close to the 
orbital plane of the comet. In such a configuration, also a neckline can become 
apparent, which is an accumulation of dust emitted 180° in true anomaly 



before the observation (IKimura and Liu 



fll977h l 



During most available observations of the 67P/Churyuni ov-Gerasimenko trail 



Ishigurd (l2008l )). the viewing geom- 



(with the exception of those described by 
etry did not favour a clear distinction between the trail on the projected orbit 
and smaller tail particles that in 3D were outside the orbit, but appeared in 
projection almost under the same position angle as the orbit. It is therefore 
difficult to discriminate between (smaller) dust emitted during the previous 



perihelion passage and (larger) dust emitted at least one orbital period before 
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the observation. From our modelling we expect that most of the dust seen in 
our images close to the projected orbit was emitted during the last perihe- 
lion passage. Nevertheless, we will in this paper use the expression "trail" for 
all dust observed close to the projected orbit of the comet, while employing 
more specific terms if referring to dust emitted during a more restricted time 
interval. 



The dust trail of comet CG has been more closely studied than most other 



comet trails, mainly beca use CG is the target o: 



European Space Agency (iGlassmeier et al. 



the Rosetta mission of the 



( 120071 )). The trail was discovered 



with the Infrared Astronomical Satellite (IRAS) in 1983 at 12, 25, and 60/im, 
and over a lengt h ranging from 0.1° in mean anomaly ahead of the nucleus 



to 1.1° behind it jSykes et al. 



(1l986al Jbh: 



Svkes and Walkeil (ll992|)). The next 



observations of the CG trail were made between Septe mber 2002 and 



ebruary 



2003 with the KISO telescope at visible wavelengths (jlshigurd (120081)). Since 
this was immediately after perihelion in August 2002, the comet was active 
and displayed a tail of young particles. The position angles of the projected tail 
and trail differed significantly, such that the old particles along the orbit were 
well separated from younger and smaller particles. In March 2003, the position 
angles of the (faint) trail of old particles and of smaller particles emitted during 
the 2002 apparition differed by a few degrees, and the old trail was not detected 
in a v isible image made A vith the Telescopio N a zionale Gali l eo on 27 March 



2003 flFuUe et al 



f l2004h : 



Moreno et al 



fl2004f )l 



FuUe et al. 



f l2004h interpret 



the observed extended dust structure in terms of a neckline and infer that the 
comet emitted a significant amount of large particles already at 3.5 AU befo re 



perihelion, which 



Kelley et al. 



las been controversial ever since ( lAgarwal et al. 



(l2007b|)). 



( 120081 ) have observed the CG trail at 3.2 AU in June 2003 (optical. 
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from Mount Palomar) and at 4.5 AU in February 2004 (24 /im with Spitzer). 
While the Palomar image was dominated by the neckline-structure, the older 
trail could be separated from the neckline in the Spitzer data. The nucleus 
of CG was closely monitored with ES0/VLT/F0RS2 and the three instru - 



ments of Spitz e r arou n d aphelion b e tween 2004 and 2007 by 



Tubiana et al. 



fl2008[ ): 



Kelley et al. 



Tubianal fcoOSf ) 



( 120091 ). Some of these observations were 



sufficiently sensitive to detect also the trail within a few arcminutes from the 
nucleus. In visible light the surface brightness of the trail was measured in 
April 20 04, June 2 004, and May 2006 at heliocentric distances of 4.7, 4.9, and 



5.6 AU f lTubianal fl2008[ ) 



Tubiana et al. 



toom . The June 2004 observation 



was performed in R, V, and I filter, resulting in the first-ever determination of 
the visible colours of the trail. In the infrared, the trail was observed at 24 /im 
in August /September 2 06 (5 .5 AU out-bound) and in May 2007 (4.8 AU in- 



bound) by 
in Table 1. 



Kelley et al. 



(120091 ). Observations of the CG trail are summarised 



[Table 1 about here.] 



In this paper we discuss data from three observations of the CG trail taken 
at heliocentric distances beyond 4.7 AU. The first observation was carried out 
in visible light with the Wide Field Imager (WFI) at the ESO/MPG 2.2m 
telescope in April 2004. It was followed in August 2005 and April 2006 by two 
observations at 24/im with Spitzer/MIPS. The comet did not appea r to be ac- 



Tubiana et al 



(l2008f ) re- 



tive during our observations, which is consistent with 
porting a point-source like nucleus during all their observations between June 
2004 and August 2006. Due to the large heliocentric and geocentric distances, 
the employed fields of view (FOVs) of about half a degree in size allowed us 
to cover a comparatively large part of the trail (up to two degrees in mean 
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anomaly). The data acquisition and processing are described in Section 2 of 
this paper. The purpose of acquiring these data was to obtain a sohd obser- 
vational basis to model the past emission of mm- to cm-sized particles from 
CG, extending the previously existing data series both in temporal and spatial 
coverage. It is problematic to infer the production rate of such large particles 
directly from observations of the coma when the comet is active, because in 
that situation large and small particles are mixed. 

Our aim was to find a model (in terms of production rate, size distribution, 
and emission speed of large particles) that would enable us to reproduce the 
observed time-evolution of the trail. To simulate images of the dust environ- 
ment for a gi ven set of parameters , we u se a generalisation of the method 



introduced by 



Finson and Probsteiru (119681 ) such that the results are valid for 



the longer timescales required to simulate a trail. The method is outlined at 
the beginning of Section 3, and the mathematical details can be found in Ap- 
pendix A. In Section 3 we also present the physical model employed for the 
cometary dust production, and introduce the variable parameters we aim to 
constrain. Simulated images for a range of parameters are shown in Section 4. 
We discuss the impact of each parameter on the morphology of the images. By 
comparison to the observations we derive the best-fitting set of parameters. 
We discuss our results in Section 5. 



2 Observations and Data Reduction 

We have observed the dust trail of comet CG at three points on the outer 
section of its orbit, once in visible light and twice at 24 ^m. Parameters of these 
observations are listed in Table 2. In this section, we describe the processing 
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applied to each of the data sets, and derive cahbrated profiles of the surface 
brightness, and the trail width along the orbit. 

[Table 2 about here.] 



2.1 Wide Field Imager Observation 

Comet CG was observed in April 2004 with the Wide Field Imager (at the 
ESO/MPG 2.2m telescope on La Silla (Chile). The heliocentric and geocentric 
distances of the comet were 4.7 AU and 3.7 AU, respectively. The camera 
has eight Charge Coupled Devices (CCDs). Each image is therefore a mosaic 
of eight frames covering a total FOV of 34'x33'. The data set includes 60 
exposures of 540 s exposure time each, of which 15 were made during each night 
between the 18 and 21 April 2004. We did not use data from the easternmost 
four CCDs from the first night because they were contaminated by straylight 
from a magnitude 4 star. The total integration time was 9h in the western 
part of the image, and 6.8 h in the eastern part. To maximise sensitivity, the 
observations were made without filter and using 3x3 on-chip pixel binning, 
resulting in a pixel size of 0.71". 



The detai 



s of 



Agarwal et al. 



fl2007ah and 



;he data processing are described in 
Agarwall (120071 ) and are briefly summarised as follows. After bias-subtraction, 
all eight exposures obtained simultaneously by the different CCDs were com- 
bined into a single mosaic, and further processing was done on these mosaics. 
Superflats were derived from the science data by normalising and median- 
combining five consecutive exposures each. Division by the superflat corrected 
also gain differences between individual CCDs and fringing effects. The fiat- 
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fielded images were airmass-corrected using a mean extinction coefficient for 
La Silla of 0.15 mag/airmass, and tlie mean sky level (approximated by the 
mode of the image) was subtracted. Object masks were obtained by average- 
combining all images of one night in the rest frame of the stars and masking all 
objects brighter than five times the local mean variance of the sky, using the 
IRAF routine obj masks. The final image of the trail was obtained by average- 
combining all images in the rest frame of the comet while applying the object 
masks and interpolating the individual exposur es to a common frame usin g 



the Spitzer Science Center's MOPEX software (IMakovoz and KhanI (120051 )). 
In the resulting image the trail is visible, but the low average signal-to-noise 
ratio precludes quantitative analysis. The image was therefore smoothed re- 
placing each pixel by the average over a rectangular neighbourhood of 140" 
parallel and 7" perpendicular to the trail. 



[Fig. 1 about here.] 



Flux calibration was achieved through aperture photometry of solar type field 
stars within the field of view of the images. Stars were considered as "so- 
l ar type" when thei r B-R and R-I filter colours in the USNO-Bl.O catalogue 



( iMonet et al. 



20031 ) were compatible with solar values. To achieve a compar- 



ison with simulated images, we ass umed a so 



sponding to a flux of 1367 W m'^ dCoxl f l2000h l 



ar R-magnitude of -27.1 corre- 



The statistical uncertainty of this calibration is 0.02 mag, which is an order 
of magnitude smaller than the uncertainty of a single entry in the catalogue, 
due to the large sample of 200 data points employed for the calibration. This 
uncertainty, however, is negligible compared to the one introduced by the 
variation in the background, which is responsible for a relative uncertainty in 
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the measured flux of up to 10% in the faint part of the trail (see errorbars in 
Fig. Ic). A systematic error was introduced by the process of superflatfielding. 
Trail information could not be completely excluded from the superflats due to 
the employed jittering pattern and because the trail signal was fainter than 
the statistical variation of the background in the raw data. To nonetheless 
compare the observation with model images, the flatfielding and smoothing 
process was included in the simulation. 



Tubianal (120081 ) has observed the CG trail within several arcminutes from 
the nucleus on 14 April 2004 with the F0RS2 instrument at ESO/VLT us- 
ing R-filter. Cross-cut profiles through the trail at 24.25" behind the nucleus 
from their VLT and our WFI observations are in good agreement around the 
projected orbit and north of it, which confirms the accuracy of the calibra- 
tions, done independently from each other. The WFI image shows additional 
brightness south of the projected orbit, which may be due to background- or 
fringing- cont aminat ion . 

The trail in our April 2004 image is characterised by a bright bulge stretching 
over 5' behind the nucleus. Beyond this region, the brightness is approximately 
constant up to the edge of the field of view nearly 35' behind the nucleus. The 
width of the trail increases steadily with distance from the nucleus, and at 
distances larger than 11' a double peaked structure is observed (Fig. Id). 
Since brightness variations of similar scales are present in the background, we 
cannot exclude that the double peak is an artefact. The difference in position 
angle of the two peaks is 0.8°± 0.2°, which agrees with the expected position 
angles of neckline and older trail (see Section 4). 

To assess whether we can detect on-going or recent cometary activity, we have 
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compared cross-cut profiles th rough the nucleus positi on with similar profiles 



through asteroids in the FOV ( iBohnhardt et al. 



(120021 )). The accuracy of this 



method is in our case low, because the pixel scale of 0.71" just barely resolves 
the seeing disc that varied between 0.5" and 2" during the observations. Within 
this limitation, we did not find any deviation of the nucleus profile from that 
of a point source and thus did not detect cometary activity. 



2.2 Spitzer/MIPS Observations 



Two mid-infrared (24 /im) observations of comet CG and its pr ojected orbit 



were carried out with the Multiband Imaging P hotometer (M I 



Rieke et al 



20041 ) on 28-29 



20041 ) on board the Spitzer Space Telescope (IWerner et al. . 
August 2005 and 8-9 April 2006, respectively. During the first observation, 
CG was at a heliocentric distance of 5.69 AU (out-bound) and at a distance 
from Spitzer of 5.72 AU. The orbit was mapped on a length of 46' projected 
to the image plane, of which 28' were in the direction trailing the nucleus. The 
observation in April 2006 took place at 5.66 AU from the Sun (in-bound) and 
5.26 AU from Spitzer. The projected length of the covered orbit section was 
42', of which 33' were in the trailing direction. 

Each observation was implemented as a cluster in photometry mode: A set of 
neighbouring pointing coordinates along the trail was specified. For each of 
these, a sequence of exposures [Data Collection Events or DCEs) was taken at 
eight scan mirror positions and repeated with a small offset of the spacecraft 
perpendicular to the direction of motion of the scan mirror. Around each 
pointing coordinate, the total field of view covered was 8.25' x 8.25', while 
a single DCE covered about 5.25' x 5.25'. The first DCE at each spacecraft 
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position had a shorter exposure time and was depressed in response. These 
images were not used for further analyses, leaving 14 valid exposures at each 
specified pointing coordinate. 

In August 2005, the set of evaluated data consists of 252 images of 30 s ex- 
posure time and each covering a field of view of 5.25' x 5.25'. Nine sets of 
pointing coordinates were specified along the comet orbit, each of which was 
covered twice, going from east to west and back again. Each section of the 
orbit - unless close to the edges of the covered area - was exposed between 28 
and 42 times, corresponding to 14 to 21 minutes in total. The images used in 
the following were processed with version 16.1.0 of the MIPS pipeline. 

In April 2006, eight pointing positions along the orbit were specified, going 
three times from east to west and back. The four easternmost pointing po- 
sitions were covered only five times. The resulting data set consists of 672 
images of 10 s exposure time and a field of view of 5.25' x 5.25'. The total 
exposure time for a given area of sky along the orbit varies between 9.8 and 
19.8 minutes, and we used the data processed with version 13.2.0 of the MIPS 
pipeline. 

From each exposure, two images were down-linked from the telescope: the 
first difference image, corresponding to an exposure of about 0.5 s, and the 
slope image, representing the flux averaged over the entire exposure time. The 
images were flatfielded and calibrated in units of MJy/sr by the automated 
pipeline at the Spitzer Science Center (SSC). A Basic Calibrated Data (BCD) 
image generated by the SSC pipeline is normally derived from the slope im- 
age. If, however, a pixel in the slope image is close to saturation, its value is 
replaced by that in the first difference image, in order to determine the fiux 
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from bright sources more reliably. The slope images from 2005 were close to 
saturation because of the long exposure time (30 s) and the bright background. 
A considerable fraction of the pixels in the BCD images had therefore been 
derived from the first difference images, such that the BCD images largely 
corresponded to images of 0.5 s exposure time instead of 30 s with accordingly 
lower signal-to-noise ratio (SNR). These images were not suitable for the anal- 
ysis of the faint trail. On the other hand, most of the pixels in the slope images 
were not hard saturated, and - as a part of the calibration pipeline - the slope 
images had been corrected for the nonlinearity of the detector response at 
high exposure levels. Therefore we used the basic calibrated slope images for 
the following evaluation. 

In 2006, the exposure time of an individual image was shorter than in 2005 
(10 s), and the background flux was lower. Hence, the slope images were not 
saturated and served as basis for the BCD images. Our analysis of the 2006 
observation is based on the BCD images. 

To better separate the cometary signal from the background interstellar medium, 
and - to a limited extent - also from stars and galaxies, shadow observations 
were carried out for both Spitzer observations, i. e. exact repetitions of the 
respective primary at a later time, after the comet had left the field of view. 

For technical reasons, the shadow for 2005 was not observed until April 2007, 
and the field of view covered by the shadow was different from that of the 
primary. No shadow is available for the trail at more than 20' distance from 
the nucleus. Since primary and shadow were made from different positions 
along the orbit of Spitzer, the zodiacal light background was different and 
had to be removed from the images before the shadow was subtracted from 
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the primary. The shadow for the 2006 observation was taken a week after the 
primary (14-15 April). 

The slope images obtained in 2005 were characterised by (1) a brightness 
gradient roughly from south-west to north-east, and (2) by a more intricate 
pattern of brightness variation that was common to all images. The ampli 



tudeof the pattern (2) diminished with increasing DCE number. 



Reach et al 



( 120071 ) remove both artefacts by subtracting a bias image from each exposure, 
obtained by median-averaging over all exposures of a given DCE number. For 
our data, different background levels in exposures taken at the same scan mir- 
ror position (same DCE number) precluded the generation of a median image 
in one step. For images of given DCE number, the background brightness in 
the primary images increases by about 0.8 MJy/sr from the easternmost to 
the westernmost images. This correlates well in amplitude and direction with 
the zodiacal light background predicted by the Spitzer Planning Observations 
Tool (SPOT). The underlying model is based on data obtained with IRAS and 
with the Cosmic Background Explorer (COBE) and takes into account the de- 
pendence of the zodiacal light background on the position of the telescope. 
The predicted contribution of the interstellar medium to the background is 
fainter by an order of magnitude. 

In both primary and shadow of the 2005 observation, we removed the zodiacal 
light and the instrumental bias by the following two steps. We first fitted a 
plane surface to the zodiacal brightness estimated by SPOT as a function of the 
pointing coordinates. We did this independently for the first and second series 
of each set of exposures to account for the time dependence of the zodiacal 
light due to the changing position of the spacecraft (about 0.1 MJy/sr over 
3 hours). The obtained plane surface was subtracted from each frame. The 
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resulting frames had approximately constant background brightness, allowing 
us to obtain their median. 

We extracted the bias image for each DCE number by median averaging frames 
of the same DCE number. For the primary, we used only frames located to the 
east of the comet, because there is no detectable trail signal in this part of the 
FOV (i.e. preceding the nucleus in its direction of motion). For the shadow we 
used all frames. The bias images were subtracted from the respective frames. 

The shadow images still showed a trend of decreasing mean brightness with 
time (about 0.3 MJy/sr over 3 hours) which would have led to noticeable arte- 
facts after averaging. Hence we fitted a linear function to the image midpoint 
as a function of the exposure sequence index and subtracted the respective 
value from each image. Thus the average background in the shadow was zero 
and its subtraction from the primary did not change the background level 
of the latter, resulting in a smooth transition between the shadowed and the 
unshadowed part (Fig. 3). 

Solar system objects are not affected by the subtraction of shadow from pri- 
mary, because they have no counterparts at the same position in the respective 
other image. They appear either white or black in the difference image, de- 
pending on whether they were present in the primary or in the shadow. Solar 
system objects were excluded from the final mosaic making use of their ap- 
parent motion between consecutive exposure series: For each observation, we 
co-added frames from the different exposure series separately (i.e. two for 2005, 
and six for 2006). By visual inspection we identified objects the position of 
which changed relative to the stars, and created mask files to flag these moving 
objects during the co-addition of all exposures. All solar system objects we 
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investigated had a non-zero apparent speed relative to the comet, thus none 
can be considered a candidate for a fragment of CG. We derive a lower size 



limit for detectable fragmen ts by comparis o n wit h the SNR of t 



le nucleus, 



which has a radius of 2 km (ITubiana et al 



l b MINI I: 



Kelley et al. 



(120091)) and 



was detected in our images at roughly 25a. If fragment and nucleus produce 
dust at the same level per surface area (if any) and have similar thermal prop- 
erties, a 5cr detection would require an object of 0.9 km radius, corresponding 



to an absolute R-magnitude of 17.1 (scalec 



the nucleus of 15.35 given by 



Tubiana et al. 



rom the absolute magnitude of 



(120081)). 



For the final mosaics we have used the zodiacal-light- and bias-subtracted 
frames from 2005, and the BCD images from the Spitzer pipeline from 2006, 
respectively. For each mosaic we average-combined all frames for the primaries 
and shadows separately in the co-moving frame of the comet in the primary, 
masking the moving objects. Each shadow was subsequently subtracted from 
the respective primary, and the results are shown in Figs. 2 and 3. The trail 
brightness and width as a function of distance from the nucleus was measured 
in perpendicular profiles. 



The surface brightness of the trail was colour-corrected to account for the slope 
of the spectrum of the dust within the spectral band of MIPS at 24 /im: To de- 
rive the true monochromatic flux density at the weighted average wavelength, 
A = 23.675 /im, the measured values were divided by a colour correction fac- 
tor of 0.948 (see MIPS Data Handbook) corresponding to colour temperatures 
between 124 K and 147 K (Section 4). The uncertainty of the colour correction 
factor due to the temperature is negligible compared with the variability of 
the background. 
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Like the optical image, both Spitzer images show a brightness bulge immedi- 
ately behind the nucleus in the sense of its orbital motion, and a comparatively 
flat profile further away. Between the three observations, the bulge moved away 
from the nucleus and became wider and shallower, both in projected distance 
and in terms of mean anomaly, indicating that the dominating particles have 
orbital periods longer than that of the nucleus. The trail width increases with 
distance from the nucleus on one flank, and is constant on the other. An in- 
terpretation in terms of different particle sizes and emission times is discussed 
at the beginning of Sec. 4. 

The images of the nucleus from both Spitzer observations have central peaks 
with a full-width at half-maximum (FWHM) of about 3 pixels (7.4" or about 
29000 km) and a surrounding diffraction ring. They are similar to other point 
sources in the same field of view, thus we do not detect any come tary activity 
at the times of our observations. However, as pointed out by 



Kelley et al 



(120091), we cannot exclude the presence of an unresolved coma within the 



central 29000 km. 



[Fig. 2 about here.] 
[Fig. 3 about here.] 



3 Simulation of Trail Images 

We outline in Section 3.1 the method used to simulate images of the comet 
trail; the detailed mathematics are given in Appendix A. Our physical model 
of the comet and the emitted dust is described in Section 3.2. 
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3.1 Simulation Method 



The position of a dust particle in a cometary trail is a function of the emission 
time and velocity, and of the radiation pressure coefficient /3, which is the 



ratio o f the solar gravity and radiation pressure on the particle flBurns et al 



( 119791 ) )■ In the following we assume that the dust emission from the nucleus is 
isotropic, that dust trajectories are determined by solar gravity and radiation 
pressure only, and that the dust properties do not change after emission from 
the nucleus. A "dust shell" describes the ensemble of locations occupied by 
dust particles characterised by a certain radiation pressure coefficient /3, and 
released at the time tc from the comet in all directions with the emission 
speed Vc- The t ime evolution of the shape and size of a dust shell is described 



as follows (e.g. 



Miiller et al. 



mm ): 



During a short interval after emission, the shell can be approximated by a 
sphere, the centre of w'hich is off-set from the comet due to radiation pressure 



( iFinson and Probsteinl . 



19681). The position r of a particle in the shell is given 



by 



r(t,/5,ve) = r(t,/3,0) + (t-te)*Ve 



where Vg is the emission velocity of the particle, and r(t, /3, 0) is the position at 
the time t of an auxiliary particle having the same /3 and tc but zero ejection 
speed. The auxiliary particle does therefore not belong to the ensemble of shell 
particles. 



On longer timescales, the shell is distorted into an ellipsoid b y tidal forces 



(IKimura and Liu 



19771 ). Equation 1 was therefore generalised by 



Fertig and Schwehm 
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fll984h to 



r(t,/3,ve) = r(t,/3,0) + $(t,te)*Ve 



(2) 



Analytica. 



Massonne 



expressions for the elements of the matrix $ were derived by 



(11985 



19871 ) from Keplerian dynamics. 



When the shell becomes more elongated, also the description by an ellipsoid 
ceases to be vahd. Large particles are then found on a tube bending around the 
comet orbit, while small particles have been dispersed by radiation pressure. 
To approximate p article positio n s on t he elongated shell, a local linearisation 



was introduced by 



Miiller et al 



(1200 ll ). For this purpose, the tube is cut into 



slices roughly perpendicular to its central axis. The auxiliary particle used 
i n Eqs. 1 and 2 is replaced by a reference part i cle in the centre of the slice 



( Kondrateva and Reznikov 



1985 



Miiller et al. 



200ll ). characterised by the 



same /3 and as the shell, but different Ve- The position of a particle on the 
slice i is expressed as 



(t,/3,ve) = r«(t) + $«(t,te,/3,v 



reft 



* Vp 



'rcf J' 



(3) 



where r^^f is the position of the reference particle at the observing time and 
v^*£ is its ejection velocity. For Keplerian motion, the elements of the matrix 



c an be calculated analytically in analogy to those given by 



Massonne 



( 11985 



19871 1 for the matrix in Eq. 2. 



We use reference particles that were emitted parallel or anti-parallel to the 
direction of the heliocentric velocity of the comet, bec ause the posit i on of a 



particle along the shell is controlled by its orbital period (jMiiller et al. 



mom . 



To first order the orbital period of the particle depends only on the component 
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of the ejection velocity that was parallel to the orbital velocity of the comet at 
the ejection time. Particles emitted parallel to the motion of the nucleus will 
eventually fall back behind it, while particles emitted backwards can overtake 
it if their emission speed was sufficiently large to compensate for radiation 
pressure. 

We choose the slices through the model tube such that they correspond to pixel 
columns in the simulated image. For every i, a reference particle (characterised 
by P and emitted at te) is determined that is located in the centre of the 
column at the observing time t. For this column, the elements of the matrix 
^^^^ are calculated. By inverting Eq. 3, the emission direction of a particle 
observed at a given point in column i is obtained. Repeating this operation 
for the four corners of a given pixel, the corresponding fraction of solid angle 
in the space of emission directions is calculated. The amount of dust located 
in this pixel follows from the overall dust production rate. By applying the 
described procedure to all pixels in all columns, the image of one particle 
shell is obtained. The final image is constructed as the sum of images of shells 
emitted at various times and characterised by various radiation pressure 
parameters /3 and emission speeds Ve- 

The trail observed along the projected orbit of 67P/C-G consists of particles 
emitted after the close encounter of 67P/C-G with Jupiter in 1959, which 
form a narrowly confined stream to date, while older particles were scattered 
into orbits very different from that of the comet (J. Vaubaillon, private com- 
munication). The dynamics of these particles relative to the comet can be 
approximated by Keplerian dynamics, while the calculation of an ephemeris 
of either comet or dust would require to take into account planetary per- 
turbations. Since we are only interested in the relative positions, we use the 
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analytical expressions for the matrix elements based on the orbital elements 
of the comet. 



The relative position of the nucleus and the dust emitted during previous 
apparitions is af fe cted by a non- g ravita tional rocket force due to outgassing 



( iMarsden 



f ll969h : 



Marsden et al. 



( Il973l )). This force affects the nucleus but 



not the dust, unless the dust was initially also outgassing. To assess the effect 
of the rocket force on the nucleus we compared ephemerides of the comet for 
the dates of our observations with and without the rocket force starting from 
osculating elements of the comet at point in time between 1959 and 2002. The 
difference between the predicted positions with and without rocket force gives 
us a measure of the error in the relative position the nucleus and dust emitted 
at the epoch of the employed osculating elements. We found that the trail 
width is affected by less than 2", which is negligible for the purpose of our 
modelling. The error of the nucleus position along the orbit is <5" compared 
to dust from 2002, <14" compared to dust from 1996, increasing up to 130" for 
dust from 1962. We expect that the distribution of dust from past apparitions 
is fairly uniform on such spatial scales (cf. Figs. 4 to 6, long-dashed lines), 
such that the shift due to the rocket force would not noticeably affect the 
model images. We therefore have not included the non-gravitational force on 
the nucleus in our model. We have here assumed that the trail particles have 



not been subject to the rocket force themselv es, i.e. have 



l ad an extremely low 



Reach et al. 



( 120091 ) have inferred 



ice content at their emission from the coma, 
this to be the case for the debris trail particles of comet 73P/Schwassmann- 
Wachmann 3 from their proximity to the nucleus' orbit. If, however, the trail 
particles had any appreciable covering of ice, it would indeed have affected 
their trajectories significantly. 
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3.2 Physical Model of the Cometary Emission 



In this section, we describe the physical model of the cometary emission on 
which our results are based. We declare, which quantities are input to our 
model, which are output, and which cannot be constrained with our method. 

In our model, dust is emitted isotropically from the nucleus. Images of the 
coma have shown az imuthal brightness v ari ations during the 1 9 96/97 and 



2002/ 03 apparitions (jSchulz et al. 



( 12004a 



Weiler et al. 



(12004) 



Schleicher 



( I2OO6I )). But no quantitative model of the distribution of activity across the 
surface of comet CG has been published to date, such that isotropic emission 
remains the simplest possible assumption. The dust emission in our model 
does not change between different apparitions of the comet after the last close 
encounter with Jupiter in 1959, because we consider the available data on gas 
production and Afp too scarce to support a more elaborate model. The comet 
in our model is active inside 3.1 AU, corresponding to a time interval of 600 
days around perihelion. We consider discrete particle classes, characterised by 
a specific value of the radiation pressure coefficient Pi, and offset from each 
other by a constant factor in (3. 

Our observations were made at 4.7 AU from the Sun and further, i.e. the 
first observation was made 10 months after the assumed cessation of activity. 
Hence, small particles with high sensitivity to radiation pressure (/3 > 0.1 in 
2004 and (3 > 0.01 in 2005/06) will have left the FOVs of our observations. For 
bulk densities of 2000 kg/m'^ and lower, all particles in the FOVs are larger 
than the wavelengths of the observed light. 

To simulate the optical observation, we assume that the light back-scattered by 
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the dust has the same spectral properties as sunhght, i.e. that the dust albedo 
is independent of (vis ible light) wavelength and particle size. In a recent paper 



Tubiana et al. 



(120081 ) have found that the trail dust has a red spectral colour, 
but since our observation was made without filter, we can only derive the 
average albedo in the visible wavelength range. The intensity incident on the 
particle is given by Jjnc = Icd/^hi with the solar flux at 1 AU Iq = 1367 W/m^, 
and the heliocentric distance of the particle rh in AU. We apply the geometric 
optics approximation, because the particles in the FOV are expected to be 
larger than visible wavelengths. The distance between particle and observer is 
A (in m). We denote the total cross-section of particles characterised by (3i in 
a given pixel by Si, i.e. Si = NiTTsj, where Ni is the number of particles and Sj 
the particle size. Our method helps to constrain Si, but neither Ni nor Si (see 
below) . The intensity received by the detector from all particles in the pixel is 

Tvis _ Pj{(^)^Q 1 r. (A\ 

■^rec - r2 A2 

where p is the geometric albedo and j(a) the phase function. 

The Spitzer images are simulated using the monochromatic flux of a blackbody 
at the weighted average wavelength of MIPS24 (23.675 /im). The intensity per 
unit frequency interval received by a detector pixel at a distance A from 
the particles having /3j is described by 

Jlf, = ^eB,iu,T)Y^S,, (5) 

i 

where B^iy^T) is Planck's function, v is the frequency, e the emissivity at 
24 jim, and T the temperature. The flux J^f^ in Eq. 5 is converted to units 
of MJy/sr by the factor lO^^/dw (pixel), where do; (pixel) is the solid angle 
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covered by the pixel. 



Since our observations are at a single infrared wavelength, we cannot infer 
the dust temperature from them. The temperature determines the intensity 
of the blackbody radiation for a given cross-section via Planck's function. 
With higher temperature, less cross-section is required in the FOV to achieve 
the same infrared flux. Assuming that the particles are characterised by the 
Bond albedo at visible wavelengths and the emissivity e at 24 /im, their 
temperature T at the heliocentric distance rh results from the equilibrium 
between absorbed solar and emitted thermal radiation: 

r(rh, Ab, 6) = 278.8 K (i^)' , ^ . (6) 



Sykes and Walkerl ( 119921 ) derived e/(l — Ab) = 0.6 ± 0.2 from IRAS measure- 
ments of the trail brightness at three wavelengths (12, 25, and 60 /im), which 
at the heliocentric distance rh=5.7 AU during our observations corresponds 
to temperatures between 124K and 147K. The IRAS-derived colour temper- 
atures are higher than the temperature of a blackbody (e=l and A-q=0) at 
the same heliocentric distance. The observed excess colour temperatures are 



consistent with randomly oriented, rapidly rotating, larg e, and dark 



capable of sustaining a latitudinal temperature gradient (ISykes et al. 



partic 



f l2004h ) 



es 



The emission speed of a particle from the coma depends to first order on 
its cross-section-to-mass ratio. We denote the inverse cross-section-to-mass 
ratio by i-s. for spherical grains of the radius s and the bul k density p. 



X = ps. Hydrodynamic models suggest that v oc (x) (e.g 



Crifo et al 
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( 119991 ): lAgarwal et al.l (l2007bf l). We use the expression 



x\ i \ 

vXo/ Up/ 



(7) 



with vq = 3.9m/s, Xo = 1 kg/m^, and the perihehon distance = 1.29 
AU. For the dependence on the hehocentric distance, a power-law with the 
exponent of -3 is assumed, because the emission speeds are to first order 
proportional to the square root of the gas production rate, and the observed 
water producti on rate can be approxi mated by a power-law with an exponent 



of about -6 (cf. 



Agarwal et al. 



(j2007bl ) and references therein) . The factor is 



a variable parameter. For = 1, the emission speeds at perihelion cor r espoii d 



et al. 



fll999h) 



to those predicted by a hydrodynamic coma model (iLandgra 
for a pe rihelion gas production rate of Qgas = 10^^ molecules/s (lAgarwal et al. 



(j2007bl )). If Eq. 7 yields a value exceeding the speed of the gas in the cited 



isotropic coma model, it is set to v 



^gas 



750 m/s. (The decrease in gas 



speed with heliocentric distance is small compared with that of the dust speeds 
and is not taken into account.) If the speed according to Eq. 7 is below the 
escape speed Vcsc, the particle is not further considered. We use Wesc=0.2 m/s. 



which corresponds to the escape speed at 20 km distance 



of 2 km (ITubiana et al. 



(l2008h : 



10^2 kg ( 


Lamv et al. 


Kellev et al. 


(2009 


)). 



( 2006 



Tom the centre of 



20071 )) and radius 



The radiation pressure coefficient (3 is characterised by the radiation pressure 
efficiency Qpr, and the inverse cross-section-to- mass ratio x- 



3 Lq Qpj. Qpi- 

P — T7 7TT7 — '^Z? ■ 

" X 



WticGMq X 

where Lq and Mq are the luminosity and mass of the Sun, and kjj = 5.77 x 
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10-^ kg/m^ fiBurns et all (11979f )l 



Both Eqs. 7 and 8 depend on x- Since v and /3 are in our model the only 
quantities determining the dynamics of a particle, each particle can only be 
characterised by its value of x-, ^lot by its size, because we have no means 
to infer the bulk density p of the dust particles. Equations 7 and 8 can be 
combined to 



/v 



Q 



--V0\ 



(9) 



pr 



The variable parameter in this equation is /v/ JQ 



;pr- 



Our method allows us to constrain the total cross-section of particles charac- 
terised by a given value of /3 emitted by the comet as a function of time. For 
constant bulk density p, the relative contribution to the total cross-section of 
particles in an i nterval (0,0 + d(3) is described by the distribution function 



f{(3) defined by 



Finson and ProbsteinI (119681 ): 



x'gix)dx oc f{(3)d(3, 



(10) 



where g{x)dx is the number of particles in the interval [x, X + dx) and (3{x) is 
given by Eq. 8. For constant p, g{x) corresponds to a differential size dis- 
tribution. We represent f{f3) by a power-law /(/3) oc (3°'. For constant p, 
the corresponding differential size distribution function is also a power-law 
dix) oc ■^^ with the exponent a = —a — A. The exponent a is a free 

parameter in our study. 

We simulate images for a range of values of f^/ ^Qpr and of a. From the simi- 
larity between simulation and observation, we determine the most appropriate 
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set of values for j^j y Qpr and a. 

We describe the cross-section production rate of particles in the interval 
d/3) for constant p by 

gs(/3,d/3,rh) = go/(/3)d/3rh-^ (11) 



Qo is a free parameter that we constrain by matching the brightness of the 
simulated trail to the corresponding observation. The dependence on heliocen- 
tric distanc e r^^ corresponds for an isotropic coma in steady state to an Afp 



parameter (lA'Hearn et al. 



f ll984l )) proportional to r^^ ^ which is in agreement 



with a power-law fit t o A^p d ata measured during several apparitions of the 



comet (lAgarwal et al. 



(j2007bl ) and references therein). 



For discrete particle classes and a given relation Sj(/3j), we have Qs,i{Pi) = 
Qi7rsf{l3i), where Qi is the number production rate of particles of size Sj. Not 
being in a position to constrain Sj(/3j), we cannot infer Qi. We can, however, 
calculate the mass production rates Qm.i of the dust for a given value of Qp^-: 



(12) 



From comparison of the cross section production rates required to reproduce 
the optical and infrared images, respectively, we infer the visible light albedo, 
because the brightness of dust at optical wavelengths is proportional to the 
albedo, while the thermal emission is to first order independent of it. Since 
the optical and the first infrared observation were separated in time by 16 
months, we can only derive the albedo on the assumption that our model 
correctly reproduces the temporal evolution of the trail. 
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We have simulated images with the original pixel scale of the corresponding ob- 
servation (0.71"/pixel for WFI and 2.45"/pixel for Spitzer). For WFI, we have 
added Gaussian noise with the parameters characterising the background of 
the observation, and simulated the fiatfielding process and the subsequent spa- 
tial filtering. The Spitzer simulations were convolved with a Gaussian having 
a = 2.3" to simulate the point spread function (PSF), and spatially averaged 
over the same ranges as the corresponding observation: 3x3 (11x3) pixels in 
the inner (outer) part of the 2005 image, and 5x5 pixels in 2006. 



4 Results 

We have simulated images for a wide range of parameters. By comparing them 
with the observations we constrain the emission speeds of dust particles as a 
function of the radiation pressure coefficient, the exponent of the dust (3 dis- 
tribution, the geometric albedo of the particles, and the dust production rates. 
Observations and simulations are compared by evaluating the peak surface in- 
tensity and the FWHM of the trail as functions of distance from the nucleus in 
plots like those shown in Figs 4 to 6. We first describe the spatial distribution 
of particles depending on their emission time. Then we assess the impact of 



separately varying the speed and (3 distribution parameters, fv/^jQpr and a. 
We derive the mass production rates and Afp parameter for the best-fitting 
set of parameters. From comparison of the optical and infrared simulation 
results, we constrain the particle albedo. 

In all three observations (Figs. 1 to 3), the distance of the half maximum 



[Table 3 about here.] 




28 



from the projected orbit is approximately constant on one side of the trail 
and increases with distance from the nucleus on the other side. The constant 
side is marked by "+" in the plots, the other by "x". The "x"-side appears in 
projection on the southern side of the trail in the 2004 and 2006 images, and 
on the northern side in 2005. Defining the "upper" side of an orbit as the one 
out of which the positive angular momentum vector points, the observer was 
below the orbital plane of the comet during the 2004 and 2006 observations 
(orbit plane angles —0.7° and —1.2°, see Table 2) and above the plane in 2005 
(orbit plane angle 1.3°). Assuming that all dust remains close to the comet 
orbital plane, the particles on the "x"-side are outside the orbit as seen from 
the inner solar system. They must therefore be subject to stronger radiation 
pressure than those remaining close to the orbit. 

In Figs. 4 to 6 the best-fitting simulations are displayed. Before describing how 
we searched for the best fit, we discuss now the spatial separation of particles 
emitted during the 2002/03 apparition of the comet on the one hand, and 
during earlier apparitions on the other. We discuss this on the example of 
the best-fit simulation, but the geometry is comparable for all other sets of 
parameters. 

Figures 4 to 6 show profiles and FWHM from four simulations: dust emitted 
at least one orbital period before the observation (1959-1997, long-dashed), 
dust emitted during the apparition immediately preceding the observation 
(2002/03, short-dashed), their sum (dotted), and the result after simulation of 
the fiatfielding for the visible image, the point spread function for the infrared 
images, and the spatial filtering for all (solid). The solid line is the profile that 
must match the observation. 
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Particles emitted 1959-1997 are concentrated around the projected orbit of 
the comet. Particles from the 2002/03 apparition form a second profile with a 
slightly different position angle, which is the cause for the broadening of the 
trail profile towards the side marked "x" in the b-panels of Figs. 1 to 3. For a 
given projected distance behind the nucleus (x-coordinate in panels a and b), 
the 2002/3 particles are more sensitive to radiation pressure than the older 
particles (Fig. 10). The difference in position angle between the two profiles 
is (-0.80 ± 0.14)° in 2004, (+1.24 ± 0.15)° in 2005, and (-0.60 ± 0.15)° in 
2006, where the uncertainty corresponds to the pixel scale employed for the 
simulations. 

[Fig. 4 about here.] 

[Fig. 5 about here.] 

[Fig. 6 about here.] 

[Table 4 about here.] 

[Table 5 about here.] 

[Table 6 about here.] 

To investigate how the shape and width of the simulated profiles depend on 
the emission speed factor J^j y^Qpr and /3 distribution exponent a, we have 
run simulations for all three observations for 105 combinations of /v/ y^Qpr in 
the range of 0.2 to 10 and a between -0.4 and 0.4. The absolute brightness in 
the simulation is controlled by the total cross-section of dust produced, which 
is a free parameter in our model. The production rates were adjusted in such 
a way that the peak brightness close to the nucleus matches the observation. 
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The employed parameter values and the results are described in Tables 4 to 6. 
The tables show a set of five entries for each parameter set. Each entry can be 
"+" (too much), "O" (acceptable), or "-" (too little). The first three of the 
five entries refer to parallel profiles like panels (a) in Figs. 4 to 6. (P) is the 
position of the peak relative to the nucleus, (F) is the brightness to the left of 
the nucleus, i.e. of dust preceding the nucleus in its orbital motion, (D) is the 
brightness in the distant trail to the right of the nucleus, trailing it. The two 
remaining entries in each field refer to the trail width as shown in panels (b) 
of Figs. 4 to 6. (T) is the width on the fiank with the old particles, (N) is the 
width on the fiank with particles from 2002/03. 

In the following, we summarise the results shown in Tables 4 to 6. We first 
describe how the choice of f^j ^JQpr and a affects the parallel profiles of the 
kind shown in Figs. 4 to 6 (a): The factor f^/ ^JQ^ mainly infiuences the profile 
of the trail close to the nucleus, i.e. the position of the peak and the shape 
of the profile in front of it. These characteristics are only little infiuenced by 
the choice of the /3 distribution exponent. When fv/\jQpi is increased, the 
bright peak close to the nucleus moves further in the direction of the orbital 
motion of the nucleus, and the brightness in front of the nucleus increases. 
Particles in front of the nucleus were emitted backwards and then overtook 
the comet, because due to their smaller total energy they are on orbits with 
smaller semimajor axis and shorter orbital period. A higher emission speed 
component anti-parallel to the direction of motion of the nucleus implies a 
larger difference in orbital period, thus the particles are found further to the 
left at a given time. Comparison of the long- and short-dashed curves in Figs. 4 
to 6 (a) reveals that most of the dust in front of the nucleus was emitted more 
than one orbital period before the observation. The sharp flank left of the 
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nucleus is due to particles from the most recent perihelion passage. The slope 
of this flank is independent of emission speed and /3 distribution, but the 
position of the peak moves to the left for increasing speed. 

Separate simulations for intervals of the considered ranges of (3 and the emis- 
sion time have shown that the point of the peak consists of particles emitted 
backwards after perihelion. All these particles have approximately the same 
distance from the nucleus. When the model speeds of all particles are in- 
creased, at the time of observation all these particles are located more to the 
direction in front of the nucleus, thus the peak shifts to the left. The particles 
in front of the nucleus were emitted around perihelion. 

The trail profile far behind the nucleus depends mainly on the /5 distribution. 
The particles found at larger distance from the nucleus are generally smaller 
than those closer to it, because they fall back more quickly due to radiation 
pressure. If we increase the relative amount of small particles by choosing a 
smaller size distribution exponent a (larger a), the relative brightness at large 
distance behind the nucleus increases. 

The width of the trail (see Figs. 4 to 6 (b)) on both sides of the orbit increases 
with the emission speed and for /3 distributions favouring smaller particles, 
the latter because smaller particles have higher speeds perpendicular to the 
orbit. 

A well-fitting set of parameters is characterised by a set of five "O" symbols 
in Tables 4 to 6 (shaded in grey). For the 2004 WFI image, good fits were 
obtained with 0.6 < fv/yfOpi < 2.2 at a = —0.2 and with 1.4 < fv/^JOpi < 
2.4 at a = 0. The Spitzer images can both be fitted with /v/ \jQpr = 1-6 and 
a = 0.2. There is no set of parameters that is ideal for both the WFI and the 
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Spitzer data. As the best compromise we adopt the values f^/ yQpr = 1.6 and 
a = 0.1. The corresponding simulations are shown in Figs. 4 to 6. 

By fitting the brightness of the simulated profiles to the observations we have 
constrained the product of geometric albedo, phase function, and cross-section 

~ 04 

production rate Qt = pj{a)Qs,i from the optical observation. The 2004 ob- 
servation was made at a phase angle of 1°, and we assume henceforth that 
j = 1. From the infrared data, we have constrained the product of emissivity, 
Planck function, and cross-section production rate = e-Bj,(i^, T)Qg j. 

We express the visible geometric albedo as a function of the dust temperature 
and emissivity: 

p = Q!''/Q''^''eB,iu,T), (13) 



which is independent of i because we assume that the albedo and emissivity 



do not depend on the particle size. Introducing Z = e/(l — A-q) and the phase 
integral q = A^/p gives p = Qi''^ /Qi'^^^^ {I —pq)ZBy{iJ, Z), where the relation 
between T and Z is given by Eq. 6. Solving for p gives the geometric albedo 
as a function of Z. In Fig. 7 we plot p{Z) derived separately from the 2005 
and 2006 observation and for two extreme values of q. We have used these 
extreme values to get a feeling for the importance of q. A commo n valu e 



Divinel f ll98lh 



would be g = 4.55, as derived from the phase function given in 
We require a geometric albedo between 0.022 and 0.044 to reproduce both 
the visible and the infr ared data with the same model parameters and in a 



manner consistent with 



were observed in dus t trails by 



comet IP/Halley by 



Ishiguro et al. 



Tokunaga et al. 



mm . 



Svkes and Walkerl (119921). S i milar ly low albedo values 



2002 



20071 ) and in the coma of 
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[Fig. 7 about here.] 

Figure 8 shows the /3-range of particles visible in the FOV of each observation 
as a function of the emission time. In general, the longer the emission lies back 
in time, the larger are the smallest particles remaining in the FOV, because of 
the action of radiation pressure. The minimum /3 (largest particles) depends 
on the assumed escape speed from the nucleus and the adopted model of 
emission speeds (Eq. 9). It is the same for all observations and apparitions. 
The maximum /3 depends on the size of the FOV, on the observation time, 
and on the assumed w(/3,rh). 

[Fig. 8 about here.] 

Figure 9 shows the minimum Afp corresponding to our best-fitting model as a 
function of the heliocentric distance. To obtain the Afp plot, we have simulated 
images for a set of observation dates taking into account for emission only 
particles within the /3-range depicted by the 2004 curve (solid) in Fig. 8. We 
used the 2004 curve because it corresponds to the largest minimum Afp of 
the three. The actual Afp at a given time may have been higher if smaller 
particles were present that have left the FOV by the time of our observations. 
The Afp-cuTves shown in Fig. 9 do not reflect so much the activity of the comet 
as the dynamics of particles moving away from the comet. For comparison, 
we plot observed Afp data: our model predicts Afp on the same order of 
magnitude or larger than the observations, which we discuss in more detail 
in Section 4. Figure 9 also shows the minimum dust mass production rate for 
geometric albedos of 0.022 and 0.044 and Qp^ = 1. The depicted production 
rates correspond to dust in the same /9-range as used for Afp. 

[Fig. 9 about here.] 
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Figure 10 shows the distribution of particle /3 coefficients along the trail for 
simulations with the best-fitting parameters. Close to the nucleus, the trail 
consists of particles having approximately 10^^ < f3 < 10"^. With increasing 
distance from the nucleus, the particles have higher /3. At a given nucleus 
distance, particles from the most recent apparition have higher /3 than from 
previous apparitions: particles in the "old trail" have /3 < 0.001, while "neck- 
line" particles are in the range of 0.001 < f3 < 0.01 in the Spitzer images, and 
0.001 < /3 < 0.1 in the 2004 WFI image. The contribution of particles with 
P > 0.01 to the neckline brightness in 2004 is 10% or less. 

[Fig. 10 about here.] 

5 Discussion 

We have investigated the properties of large dust particles emitted by comet 
67P/Churyumov-Gerasimenko through analysis of three images of its dust 
trail - one in visible light and two at 24 /im - obtained at heliocentric dis- 
tances between 4.7 AU and 5.7 AU. We have generated simulated images of 
the cometary trail, and derived model parameters suitable to reproduce the 
observations. 

Assuming that the comet is not active outside 3.1 AU, the dust present in 
the images is characterised by radiation pressure parameters /3 < 0.01 in 
the infrared images, and /3 < 0.1 in the visible image. For a bulk density 
of 1000 kg/m^, this corresponds to minimum sizes of 60 fim and 6 fim, respec- 
tively. Our data do not allow us to make any conclusions about grains having 
higher (3. 
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Our model has four free parameters: the relation between emission speeds 
and radiation pressure coefficient /3, the exponent of the dust /3 distribution 
(or - for size- independent bulk density - differential size distribution), the 
production rate of dust cross-section, and the average geometric albedo of the 
dust at visible wavelengths. 

The closest reproduction of the observations was achieved with a /3 distribution 
exponent of a = 0.1 (differential size distribution exponent a = —4.1), and 
fv/ \jQpr = 1-6 in Eq. 9. A particle having /3 = 0.01 is emitted with 25m/s 
at perihelion and with 2m/s at 3 AU. The production rates of particles larger 
than 10 /im are on the order of 1000 - 2000 kg/s at perihelion and of 1 - 2 kg/s 
at 3AU. 



Our model repr oduces the su r: 



February 2004 (IKellev et al 



ace b rightness of the trail seen with Spitzer in 



( 20081)') and of the trail profiles obtained with 



VLT in April and June 2004 (ITubianal (120081 )) within a factor of two. It also 



reproduces the surface 



mm; 



Moreno et al. 



jrigh tness of coma and tail in March 2003 (iFuUe et al 



(120041 ) )■ if taking into account only particles in the /3 
range depicted in Fig 8. This implies that our model does not allow for ad- 
ditional, smaller particles being emitted from the comet, which at that time 
was active. We h ave simulated our images with the parameters derived by 



Fulle et al 



(j2004[ l (see Fig. 11) and found that their model underestimates 



the surface brightness in our data by an order of magnitude. 



The production rates, emission speeds, and /3 distribution derived in Section 4 
correspond to Afp on the order of 300 cm at perihelion and 50 cm at 3 AU post- 
perihelion. These values include only particles in the /3 range shown in Fig. 8, 
i.e. having /3 < 0.1. Thus, with our model, the production of large particles 



36 



alone leads to Afp on the same order of magnitude or larger than the measured 
Afp. The P distribution exponent of a = 0.1 implies that in our model the 
largest contribution to Afp stems from the smallest particles present. Since 
the particles with /3 < 0.1 more than reproduce the observed Afp, we conclude 
that the exponent a = 0.1 cannot be valid for smaller particles. Particles with 
P > 0.01 are in our data only present at distances larger that 5' behind the 
nucleus in the 2004 image. In our best-fitting model they contribute about 10% 
to the flux from this region (Fig. 10), and the model slightly overestimates the 
total flux (Fig. 4). If the emission of particles with 0.01 < /3 < 0.1 was less 
than given by the assumed power-law, we expect that both the simulated Afp 
in Fig. 9 and the simulated trail in Fig. 4 would be closer to the observations, 
while the simulations of the Spitzer images would not change. We conclude 
that the coma brightness while the comet was in the inner solar system may 
have been dominated by particles with /3 ^ 0.01 (60 - 600 /im, for densities 
between 1000 and lOOkg/m^). 



Hanner et al. 



( 1l985l ) in thermal infrared spectra of the CG coma have found 



excess colour temperatu res that are consis t ent ^y ith the temperature excess 



observed in the trail by 



descri bed by Eq. 6 and Z = 0.6 ± 0.2 (lAgarwal et al 



Sykes and Walker! fll992[). AH measu re ments can be 



(120091)) 



Divine et al 



( 1l986l ) explain the excess colour temperature in the coma with a dominant par- 
ticle size of 0.6 pm. Alternative explanations could include grains sufficiently 
large to sustain a temperature gradient, or a certain degree of porosity. 



A different explanation for the apparently high production rates of large par- 
ticles required to reproduce the trail images, would be the fragmentation of 
large particles several years after their emission. If the trail particles had still 
been disintegrating at the time of our observations, there would have been a 
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population of small particles in the FOVs that may have contributed signifi- 
cantly to the measured surface brightness. Thus the amount of large particles 
originally produced from the comet could have been smaller, leaving room for 
a contribution of smaller particles to the coma flux. The cause of a possible 
fragmentation remains at present speculative. 

Although we have assumed isotropic dust emission from the comet, our method 
primarily constrains the component of the emission velocity perpendicular to 
the Sun-comet axis: The position of a grain relative to the nucleus along the 
orbit depends mainly on its velocity component parallel to the orbital motion 
of the nucleus at the emission time, while its distance from the projected orbit 
depends on the velocity component perpendicular to the comet's orbital plane 
and on the time elapsed since emission. The dust emission speeds close to the 
subsolar point are therefore less constrained by our model. We expect that an 
insolation-driven emission model, where the local production rates and speeds 
depend on the solar zenith angle, will also reproduce the observations, as long 
as the velocity components perpendicular to the comet-Sun vector are similar 
to those we found with the isotropic model. A different possible non-isotropic 
scenario could include high emission rates of particles with a large velocity 
component parallel to the direction of motion of the comet. Such particles 
would eventually fall back behind the nucleus and could form the observed 
bulge. 

To put our findings in the context of published models we show in Fig. 11 a 
comparison of the results from different models for the dust size distribution, 
P range, emission speed, and dust production. Details of the observations 
on which the different models were based are listed in Table 1. The time- 
dependence of all parameters is a free parameter for 



Pulleetal 



(120041 ) and 
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Moreno et al. 



(120041 ) ■ while 



Ishigurol (I2008f l: 



Kelley et al. 



( 20081 



20091 1 and our 



model assume variable but time-independent size distributions and power- 
laws with heliocentri c distance wit h var iable exponents for emis sion speed 



and dust production. 



Ishigurol (120081) and 



distribution exponent of -3.5, for 



Ke 



FuUe et al 



lev et al. 



(12008 



20091 ) find a size 



( 120041 ) it drops from -3.5 before 



perihelion to -4.4 after. Their intermediate values around perihelion are s imilar 
to our result of -4.1. In all models - with the exception of 



Moreno et al. 



mom 



the dust cross-secti on is domi n ated b y the smallest parti cles (a < —3). 



The P range given by 



FuUe et al. 



(120041 ): 



Moreno et al. 



hoo4} and us refiects 



mainly the dates and FOVs of the modelled observations, with an additional 
infiuence of the emission speed (particles emitted at higher speed parallel to 
the orbital motion disappear more quickly from a given FOV). In our model. 



the largest grain size 
and the adopted v{(3). 



depends on th e ass umed escape 



Ishigurol ( 120081 ) and 



Kelley et al. 



speed of t 



le comet 



(12008 



2009|) seem to 



assume fixed size ranges: due to the shallow dependence of the emission speed 
on heliocentric distance in their models {r^^'^), the speed of a particle with 
/3 = 10~^ is larger than the es cape speed even at ap helion. The emission speeds 



in all models except that of 



Moreno et al 



( 120041 ) are rough! 



peri helion, but differ by about an order of magnitude at 3 AU. 



and 



y consistent at 



Kelley et al. 



(12008 



llshigurol (I2OO8 



20091 ) assume v oc r^^'^, which represents the varying 
insolation with rh. We use v oc r ^^ based on the observed gas production 



rate of the comet. In the model of 



FuUe et al 



(120041 ) the time dependence of 



the emission speed is a free parameter, and they find a strong decrease from 
before to after perihelion, which would indicate that also the gas production 
should be asymmetric about p erihelion. Note that v oc /3^^'^ in all models 



except for 



Moreno et al. 



( I2OO4I ) who use a variable exponent and find fitting 



solutions also with v oc l3^^^. Similar to the speeds, the dust production rates 
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are free to vary with time in 



Fulle et al. 



find a strong asymmetry with perihehon. 



fl2004h and 



2009) 



Ishigurol ( 12008 



Moreno et al 



(200 



Kelley et al 



) , who 



( 12008 



ana 



our m odel assume a power-law dependence on heliocentric distance. 



Ishigurol ( 120081 ) tre ats the exponent a s a va riable parameter and finds -3 to 



fit his images best. 



Kelley et al 



(12008 



20091 ) assume an exponent of -5.8, but 



do not give absolute production rates. Our model assumes an exponent of -8. 
The dust product i on pe r orbit in our model is about 2.5 times higher than 



found by 



Ishigurol (120081). 



[Fig. 11 about here.] 



We wish to stress that in this paper we have presented one possible model that 
reproduces the discussed three observations and does not entirely fail for a set 
of others. We have assumed fixed dependencies of emission speed and pro- 
duction rate on the heliocentric distance, a constant /3 distribution exponent, 
isotropic emission, and size-independent optical properties. We have done this 
because the observational evidence to support more elaborate models is scarce. 
However, models with different assumptions on e.g. the dependencies on he- 
liocentric distance and the e mission geo r netry r eproduce the respect i ve obser- 



vations equally successfully 



(120081); 



Kelley et al. 



( 12008 



Ful 



e et al. 



(|2004h : 



Moreno et al 



( l2004l ): 



Ishiguro 



20091 )). It must be the task for future work to find 
ways of discriminating between the differences in the existing models. A key 
to such effort lies in the increasing amount of observations of comet CG, char- 
acterising both the time-dependent activity of the nucleus and the temporal 
evolution of the dust trail. 
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A Mathematical method 



We describe in this appendix the algorithm to simulate the image of a shell 
of particles characterised by the radiation pressure parameter /3, the emission 
time tg, and the isotropic emission speed Vg for a given observation time and 
observer location. 



A.l Coordinate Systems 



The coordinate system used to describe the image has its origin at the comet 
nucleus, and {x, y, z) form a right-handed, orthogonal system. The 2;-axis is 
along the line of sight, pointing away from the observer. The x-axis is parallel 
to the Earth's equatorial plane, and x increases with right ascension (RA), 
from right (west) to left (east) for an observer looking at the sky. The y- 
axis completes the orthogonal system, pointing northward in projection. It is 
assumed that the x-y-plane corresponds to the image plane, which introduces 
a distortion for larger fields of view. 



A second co met-centred coo rdinate system is the Cometocentric Bipolar Sys- 



tem (CBS) f lMassonnd . 



19871 ). It is spanned by the orthogonal vectors (^,X)C)) 
where the (^-axis is along the line Sun-comet, pointing away from the Sun. The 
^-axis is in the orbital plane of the comet, perpendicular to C, and pointing in 



the direction of motion, x completes the right-handed system. 
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A. 2 Reference Trajectories 



A set of reference particles is defined that have the same /3 and as the shell 
to be depicted, but not the same emission speed. All reference particles have 
been emitted in the direction of the heliocentric velocity vector of the comet 
Gv.nucK^c), with varying absolute speeds. One reference particle J is needed for 
each pixel-column Xj in the image. To find the corresponding emission speeds 
v\\j{xj), an auxiliary set of particles i with emission speeds v\\^i is defined and 
the particle positions Xj in the image calculated. The speeds of the reference 
particles v\ij{xj) are found by interpolation of the relation Xi{v\^,i). The orbital 
elements of each reference particle j are calculated from its heliocentric state 
vector at emission and /3. The coordinates {xj, Uj, Zj) of the reference particles 
at the time of observation are denoted by the vectors roj. 



A. 3 Linearisation 



The position r(tobs) of a particle with emission velocity Vc = i'||,jev,nuci + Av 
is linearised around roj for not too large Av: 

r(tobs) = roj(tobs) + *jAv. (1) 



The transformation matrix depends on the trajectory of the reference par- 
ticle and on the dates of emission and observation, ^'j is defined such that Av 
is given in the CBS-frame, whereas the positions r are in image coordinates. 



If the particle trajectories are considered as Keplerian orbits with the gravi- 
tational potential modified by the factor A/x = (1 — th e elements of the 



matrix are obtained in analogy to iMassonnd (119871 ). While 



Massonnel fllQSTI ) 
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uses the nucleus as reference object, we use d ust particles. The full transfor- 
mation formula is as follows ( iMassonnd . Il987l . pp. 127-128): 



.CBS 



it) 



6T/5t\t 



( \ 

$22 $21 
$33 



^$12 $iy 



afl-e' 



(2) 



Part II describes the change in trajectory due to the different strength of 
radiation pressure for the nucleus and a dust particle. In the present approach, 
/3 is the same for the reference particle and the shell particles. Hence part II 
is always zero. In part I, bTjbt\t^ is the time derivative of the true anomaly of 
the reference object at the time of emission: 



\ a3(l-e2)3 



1^1 + e cos Tc)^ 



(3) 



/i = GMq = 1.33 X 10^° m^/s^ is the gravitational potential of the Sun, and 
a, e, and r are the semi-major axis, eccentricity and true anomaly of the 
reference object, respectively. /3 is the radiation pressure coefficient of the 
reference object (i.e. /3 = in the situation described by Massonne). The 
matrix elements $mn are functions of tobs, ^e, Tohs, To, a, and e. r*"^^(t) and 
Vjp^^(tc) are in CBS-coordinates centred on the reference object at the times 
t and te, respectively. 
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The Av of Equation 1 is equivalent to v^, (tc) of Equation 2. It follows that 



Tcbs— im/j \ 
j V'ohs J 



$22 $21 
$33 



$12 $11 



(4) 



where Tj^''^~™(tobs) is the matrix to transform reference-particle centred CBS- 
coordinates to image coordinates relative to the reference particle. The orbital 
elements in ^'j and 5Tj/5t\t^ are those of the reference particle, r,-, Oj, and Cj. 



A. 4 Reference Frame in Ejection Velocity Space 



A reference frame in the space of ejection velocities is defined with unit vectors 
Uj, Vj, Wj such that 



fie. 



UjVy = 



(5) 
(6) 
(7) 



(9) 



A particle emitted with Av = uj will be located on the same line of sight as 
the reference particle, i.e. a; = Xqj and y = y^j. For Av = v^, the particle 
will have x = xqj, but y ^ yoj. To be consistent with the definition of ^j, the 
unit vectors uj, Vj, Wj need to be given in the CBS-frame at the time tc. By 
Av = puj + aVj, a plane in the Av-space is defined that contains all possible 
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Av for particles to lie in the image column characterised by x = xoj. If 



uvw— cbs 



„,cbs „,cbs „,,cbs 



„ cbs cbs „,,cbs 
"j,2 ^j,2 "^j,2 



cbs cbs „,,cbs 
"j,3 ^j,3 Wj,3 



(10) 



is the matrix to transform coordinates from the (M,w,w)-system to the CBS 



frame, the matrix ^\ = ^\Tl 



uvw— cbs 



^^^^ iiiu,^xi^v — ^j-i-j by definition takes the form 



^ ^ rji uvw— cbs 

J J J 



ci 

Z/ C2 



where the vector c = ^'jWj has been defined, c, A, fi, and u depend on j, but 
the index is omitted here for the sake of legibility. 

A. 5 Mapping Emission Direction to Position in Image 

Any vector v can be expressed as a sum of the unit vectors uj , \j , Wj : 



v = V sin 6 sin u,- + f sin 6 cos v ,■ + w cos 6 w , 



(12) 



where ^ G [0, vr] and G [0, 27r]. Hence, Av is 

Av = Vc - t'llJ ev,nucl 

= Ve sin 9 sin Uj + Vc sin 6^ cos Vj + cos 9 wj — jev,nuci- (13) 
Taking into account Equations 5-9, Equation 4 reads 
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r = ro,j + *jAv (14) 
= roj- + Vc {sm9 sin0 *jUj+ sm9 cos0 ^jv^ + cos^ ^jWj) — f y j ^je^ nuci 

Ae^ 1/ey+iJ.ez c cr 

or component-wise 



x = xoj + VeCicosO — v\\jai (15) 
y — yo,j + "^c C2 COS 6^ — t>|| J + f c sin 6^ cos 0. (16) 

(9) 

These equations relate the direction of emission of a particle to its position in 
the image. They can be inverted to give 6{x) and (f){x,y): 



cose='^M^l±^^^^ (17) 

VeCi 



y-yj{e{x)) 

cos = . 

Ay^ax(^(x)) 



For fixed 6*, the full ring-segment in a pixel-column covers all values of G 
[— 7r,7r]. The borders of the projected shell are characterised by = vr and 
= 0, respectively. Hence, the two parts of a shell (front and back) in the 
given pixel are described by the two angles 0(^) and —(f){9). The sign of sin0 
does not change within a given side of the shell. 



A. 6 Parameters of the Map 

To exploit Equations 17 and 18, the parameters u, ci, C2, ui, and (T2 need to 
be evaluated. First, the inverse matrix is calculated. The column vectors 
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in ^'j ^ are denoted as ^j(i) with i = 1 . . . 3 such that 



*;' = (*7d)'*7(2)'*7(3)) 



*7i!(i) *7i!{2) *7i!(3) 



*79^1^ ^79(9\ *79^^^ 

j2,(l) j2,(2) j2,(3) 



y ^i3,(l) ^j3,{2) ^i3,(3) J 



(19) 



Equation 5 yields that uj = X^j^^y With |uj| = 1 (Equation 7) follows 



:*7i!(3))^ 



;*72!(3))^ + (^73!(3))^ 



(20) 



Equation 6 gives 



V, = ^*-^)+/^*7(3)- 



(21) 



Multiplication with Uj results in 



= U, V, = A ^^.(l^) *i(3) + ^ (*J(3))'> 

1/A2 



(22) 



which can be solved for /i/z/ (A, ^j"^)- Re-insertion into Equation 21 gives the 
vector "Vj/u, from which u is obtained due to the normalisation of Vj. Then it 
is straightforward to calculate Wj = Uj x v^, c = ^'jWj, and cr = ^'j nuci- 
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A.l General Solution 



The aim of this section is to calculate the number density of particles in a 
given pixel confined by the coordinates (xi, yi) and (x2,?/2)- The distribution 
of particles released to a unit solid angle is given by 



dn^'^'^^-^d^' ^^^^ 



where 6{x) and (f){x,y) are the comet-centred angles defined in Equations 17 
and 18. The amount of dust in this pixel is obtained by 



X2 y2 



f f dN 
N= / / ——dydx 
J J dxdy 



dy 



x\ yi 

00362(2:2) 4>2iy2,cosd) 
dcos 6 



cos9i(a:i) (^i (j/i ,cos ^?) 



dN 
dlT* 



, cos( 



+ 



00862(2:2) — <^i(j/2,cos6) 



f dN 
dcos6' J d(p -^{(l),cos9) 



cos6i(a;i) —<f>2{yi,cosl 



|iVl| + |iV2|. 



(24) 



The first integral represents the number density in the side of the shell with 
< < vr while the second integral represents the side with negative (p. 
Therefore, Ni and N2 cannot compensate each other. 

If a horizontal boundary of the pixel is inside the projected shell, the corre- 
sponding limit of the inner integrals in Equation 24 is given by Equation 18. 
This is in the following referred to as a "normal" integration limit. If the edge 
of the projected shell is inside the pixel, the limit ±0j {yj,cos6) of the inner 
integral must be set to = or = tt as appropriate (henceforth called an 
"anomalous" integration limit). Introducing z = cos 6, the inner integration 



48 



limits in Equation 24 are described by 

r 

<Pnmiy,z) I COs(0nm)| < 1 

(l)(y,z) = < 

(h^r, else 



(25) 



with (f>nni{y, z) given by Equation 18, and </)an G {0, tt}. The careful evaluation 
of contributions from the edges of the projected shell is critical, because a 
significant fraction of the dust can be concentrated in these pixels due to the 
shallow angle between the line of sight and the surface of the shell. 

[Fig. 12 about here.] 

As an example, the integral A^i is given for the situation depicted in Figure 12. 
Introducing dN/dfl = f{4), z) and its antiderivative F(0, z), the integral A''i is 

4>2{y2,z) 



m = Jdz J d<pf{<p,z) 



1 <l>i{yi,z) 

22 

dz [F{4>2{y2,z),z)-F{Myi,z),z)] 



(2) 



dzF{(j)nmiy2,z),z) + / dzF {(f)^n{y2, Z2) , Z 



^0 



^0 



-J dzF{(j)^niyi,Zi),z) - J dzF{(f)nmiyi,z),z). 



(26) 



^0 



For each horizontal pixel boundary, five scenarios of intersection with the edge 
of the shell can be distinguished. These are depicted in Figure 13. 



[Fig. 13 about here.] 
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A. 8 Isotropic Dust Emission 



Isotropic emission is described by / = dN/dQ = Q/An. In this instance, the 
integral in Equation 24 can be solved analytically. For symmetry reasons, the 
contributions from the front and back sides of the shell are identical, and 



Nt"^ = Jdz /d0 ^(0, z) = ^ldz [0(2/2, z) - z)] = Nt^ - Nt\ (27) 



where Zi = z{xi), = z{x2), and "a" and "b" refer to the upper and lower 
pixel boundary {y = 1/2 and y = yi), respectively. 

Normal integration limits. If the part of the upper pixel boundary that is 
inside the shell is limited by z^ and Zjs (in the example given in Equation 26 
Za = z\ and z^ = z^''), the corresponding normal part of the integral A^'^^"^ is 
given by 



Equation 24 simplifies to A^^'*^") = 2 ivf with 



Z2 0(?/2,z) 



21 </>(yi,2) 



N 



(iso) 



Q 



I 



dz 0nm(2/2, Z) 



la 
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with 



1 



H — arctan , , , , 

2 V {A2 + B)^2 



A, 



arctan 



P2,2 



(28) 



A2- 

B 

Pl,2 ■■ 
P2,2- 



2/2 - Vjfi + vija2 



02_ 
V 



{A2 + Bz) 



z + B'z + A2B. 



(29) 

(30) 

(31) 
(32) 



The second index "2" indicates that A2, pi,2, and ^2,2 depend on the upper 
pixel boundary {y = 1/2). 

If the edge of the shell intersects the pixel boundary, one or both of Za and Zjs 
in Equation 28 are given by an intersection point rather than by a pixel corner. 
As discussed at the end of Section A. 5, for the value of z in question follows 
cos 01/2 = ± 1, and some simplifications can be introduced into Equation 28. 
To avoid numerical errors at the intersection points, it is advisable to introduce 
these simplifications explicitely into the computer code. 



A2 + BZa/l3 



±1 



-a/13 
Pl,2 = 0, 

P2,2 = Za/l3 + hosBsine, 



arctani 



G 



/Pia 



sign(G). 



(33) 

(34) 
(35) 

(36) 



The expression G/^/pi^ refers to any of the three arguments of the arctan in 
the last three lines of Equation 28. The simplified expression for the concerned 
term in N[^a°a/i3 then reads 
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/\r(iso) _ Q 



kz, 



«//3 



TT 

- ^ Sign 



- B) 



TT . ( Al + A2Bz^/p-l+P2 

H — sign ; 



- sign (P2) 



(37) 



The solution of n'^^°^ is analogous to Ni^°\ replacing ?/2 by yi, and in general 
with different values for and 2;^. In the example given in Equation 26 the 



lower boundary would be characterised by z^ = Zq'^ and Zjs 



Anomalous integration limits. For those sections of a pixel boundary that 
are outside the shell (limited by z^ and zs), the inner integration limit in 
Equation 27 simplifies to (f>{yi, z) = </)an = ki, where the index "i" refers to the 

upper and lower pixel boundary (?/2 or yi), and the value of ki is either or 

(2) (3) 

IT. In the example given in Equation 26, z^ = Zq ' and zs = Zq for the upper, 
and = z\ and = z^"^ for the lower pixel boundary. The anomalous part 
of the integral N^^^jl Equation 27 reads then 



Nt^ = ^ Jdz My2,z) 

Q 



As for normal boundaries, the solution for the lower pixel boundary, N[^^°^ , is 
analogous to Nia°\ replacing y2 by yi, and having different values for z^ and 

Zs. 
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Table 1 

Observations of the CG trail and/or neckline. The heliocentric distance is given 
in AU. "Kiso" refers to the Kiso Observatory 1.05m Schmidt telescope at Nagano, 
Japan; "TNG" is the Telescopio Nazionale Galileo at the Observatorio Roque de los 
Muchachos on La Palma, Spain; "Palomar" refers to the Mount Palomar 5m Hale 
telescope Large Format Camera; "MIPS" is the Multiband Imaging Photometer 
for Spitzer; "F0RS2" is the FOcal Reducer and low dispersion Spectrograph at 
ESO/VLT; and "WFI" is the Wide Field Imager at the ESO/MPG 2.2m telescope 
on La Silla. AMA is the range in mean anomaly corresponding to the fields of view 
of the observations, except for IRAS, where the range of the detected trail is given. 



Date 



Th Instrument Wavelength AMA [°] Reference 



May 1983 2.3 IRAS 



12,25,60 /xm 0.1 .. 



Sep 2002 


1.3 


Kiso 


R 


0.0 


Dec 2002 


1.8 


Kiso 


R 


0.1 


Feb 2003 


2.2 


Kiso 


R 


0.1 


Mar 2003 


2.6 


TNG 


R 


0.0 


Jun 2003 


3.2 


Palomar 


R 


0.3 


Feb 2004 


4.5 


MIPS 


24 /xm 


0.1 


Apr 2004 


4.7 


F0RS2 


R 


0.0 


Apr 2004 


4.7 


WFI 


vis, no filter 


0.1 


Jun 2004 


4.9 


F0RS2 


R,V,I 


0.1 


Aug 2005 


5.7 


MIPS 


24 fxm 


1.0 


Apr 2006 


5.7 


MIPS 


24 /xm 


0.5 


May 2006 


5.6 


F0RS2 


R 


0.1 


Sep 2006 


5.5 


MIPS 


24 yum 


0.2 


May 2007 


4.8 


MIPS 


24 fxm 


0.1 



-1.1 Svkes et al. (1986b^), 

Svkes and Walker (1992) 
-0.2 Ishiguro (2008) 
-0.3 Ishiguro (2008) 
-0.2 Ishiguro (2008) 
-0.1 Fulle et al. (2004), 

Moreno et al. (2004) 
-0.6 Kellev et al. (2008) 
-0.4 Kellev et al. (2008 ^ 
-0.1 Tubiana (2008), 

Tubiana et al. (2008) 
-1.0 Agarwal et al. (2007a), 

this paper 
-0.1 Tubiana (2008), 

Tubiana et al. (2008) 
-1.6 this paper 
-1.7 this paper 
-0.1 Tubiana (2008), 

Tubiana et al. (2008 ^1 
-0.2 Kellev et al. (2009) 
-0.1 Kellev et al. (2009) 
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Table 2 

Geometrical parameters of the three observations of the dust trail of comet CG. 
MA = mean anomaly, PA = position angle. The "orbit plane angle" is the an- 
gle between observer and comet orbital plane, measured from the centre of the 
comet nucleus; positive values indicate that the observer is above the comet or- 
bital plane, in the direction of the orbital angular momentum vector of the comet. 
The "viewing angle" is the angle between the line of sight and the comet veloc- 
ity vector; values < 90° indicate that the comet moves away from the observer. 
The "neckline emission date" is defined as the date when the comet nucleus was 
at a true anomaly of Tobs ~ 180°. With the exception of the orbit section lengths 
and the viewing angle, all values were obtained from the JPL Horizons System 
(http : //ssd . jpl . nasa . gov/horizons . cgi). 



Telescope 


ESO/MPG 2.2m 


Spitzer 


Spitzer 


Instrument 


WFI 


MIPS24 


MIPS24 


Date of observation^ 


18-21 Apr 2004 


28/29 Aug 2005 


8/9 Apr 2006 


ITelinceTitrir distanre 


4.69 AU 


5 69 AU 


5 66 AU 


Distance from observer 


3.69 AU 


5.72 AU 


5.36 AU 


Orbit section covered (proj.) 


2'. . . -33' 


18'. . . -28' 


9'. . . -33' 


Orbit section covered (MA) 


0.1°... -1.0° 


1.0°... -1.6° 


0.5°... -1.7° 


PA of neg. velocity vector 


296.9° 


293.5° 


281.6° 


PA of Sun-comet vector 


253.1° 


106.2° 


274.9° 


Phase angle 


1.0° 


10.3° 


10.0° 


Orbit plane angle 


-0.7° 


1.3° 


-1.2° 


Viewing angle 


54° 


94° 


110° 


True anomaly at obs. (robs) 


150.8° 


176.1° 


185.9° 


Neckline emission date (tni ) 


15 Jul 2002 


14 Aug 2002 


25 Aug 2002 


tni relative to perihelion 


-34 d 


-4 d 


+7d 


Heliocentric distance at t^i 


1.36 AU 


1.29 AU 


1.29 AU 



^The values given in this table are for 20 Apr 2004, 29 Aug 2005, and 9 Apr 2006, 
00:00:00.0 UT. 
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Table 3 

Free parameters of the employed model. Columns 1-3 summarise the physical 
model outlined in Section 3, column 4 describes the main effect (s) a parameters has 
on the simulated images (see Section 4). 

Physical quantity Model Variable Related observational quantity 

Diff. /3 distribution Eq. 10 a Brightness far behind nucleus 

relative to peak; Trail width. 

Emission speed Eq. 9 f^/ ^JQp^ Profile close to nucleus 

(peak position and profile 
in front of peak); Trail width. 

Cross-sec. prod, rate Eq. 11 Qq Absolute brightness of trail 

in IR for given temperature. 

Geometric albedo P P i^-, s) p Relative brightness in IR and visible. 
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Table 4 

Simulation results for the WFI observation in April 2004. We have run simulations 
for 21 values of the parameter /^/ \/Q^ and five values of the /3-distribution ex- 
ponent a. For each pair of parameters we have plotted profiles like Fig. 4 and by 
visual inspection compared the simulation and observation with respect to five cri- 
teria: "P", "F", and "D" were evaluated in a profile like in panel (a), and "T" and 
"N" refer to the trail width as in panel (b). P is the position of the peak bright- 
ness near the nucleus. "-" ("+") means that the peak is too far behind (in front 
of) the nucleus in the direction of its orbital motion. F characterises the brightness 
in front of the nucleus (positive x- values). "+" ("-") indicate that the trail in this 
region is too bright (faint). D evaluates the trail brightness in the nearly flat part 
of the profile far behind the nucleus. "+" ("-") have the same meaning as for "F". 
T and N refer to the width on the trail ("+") and neckline ("x") side of panel 
(b). "+" ("-") means that the simulation is wider (narrower) than the observation. 
In all cases, "O" means that the simulation reproduces the observation within the 
errorbars. "n.a." for P indicates that the profile is too flat to identify a peak. Good 
fits are shadowed in grey. 
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Table 5 

Simulation results for the Spitzer observation of August 2005. For details see Table 4. 
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Table 6 

Simulation results for the Spitzer observation of April 2006. For details see Table 4. 
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Fig. 1. Optical image of comet CG obtained in April 2004 with the Wide Field 
Imager at the ESO/MPG 2.2m telescope on La Silla. The four panels are aligned in 
X, the scale in projected distance from the nucleus is given at the bottom of panel 
d, and the corresponding offset in mean anomaly of the comet is shown at the top 
of panel c. One arcminute corresponds to l.GxlO^m at the comet, a) Unfiltered 
mosaic image, b) Same image, after replacing each pixel by the average over a 
neighbourhood of 200 pixels (140") parallel and 10 pixels (7") perpendicular to 
the trail axis (boxcar-filtering). The size of the filtering window is indicated in 
the upper right corner. The data points in panels (c) and (d) were derived from 
profiles perpendicular to the projected orbit at discrete distances x from the nucleus. 
Examples of such profiles are displayed in panels (c) to (f) in Fig. 4. c) Height 
above background of the brightest point in the trail cross-section. The x-errorbars 
are the size of the filter window. The y-errorbars are the rms of the background. 
Dashed lines indicate the corresponding surface brightness in R-filter. d) "+" and 
"x" show the distances from the orbit where the surface brightness is half of the 
peak brightness, "+" indicating the trail fiank that is closer to Sun and observer. 
Circles show the peak position when only one peak is observed. Otherwise, squares 
and diamonds indicate the positions of the peaks attributed to old trail and neckline, 
respectively (see beginning of Sec. 4). 
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Fig. 2. Spitzer/MIPS24 observation of the CG trail on 28-29 August 2005. All 
panels are aligned in x, the scale in projected distance from the nucleus is given at 
the bottom of panel d, and the corresponding offset in mean anomaly of the comet 
is shown at the top of panel c. One arcminute corresponds to about 2.5xl0^m 
at the comet. One pixel corresponds to 2.45". a) Mosaic of the primary images. 
The surface brightness is in MJy/sr, and moving objects close to the comet orbit 
have been removed, b) The same image after subtraction of the shadow taken in 
April 2007. To derive the data points in panels (c) and (d), we boxcar-filtered the 
image in panel (b) over 3x3 pixels close to the nucleus and 11x3 further away, and 
evaluated profiles perpendicular to the projected orbit at discrete distances x from 
the nucleus. Examples of such profiles are displayed in panels (c) to (f) in Figs. 5. 
c) Height above background of the brightest point in the trail cross-section. The 
error bars in x are the size of the boxcar averaging window, those in y represent la 
of the background variation, d) FWHM of the trail (see caption of Fig. 1). 
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Fig. 3. Spitzer/MIPS24 observation of the CG trail on 8-9 April 2006. All panels 
are aligned in x, the scale in projected distance from the nucleus is given at the 
bottom of panel d, and the corresponding offset in mean anomaly of the comet is 
shown at the top of panel c. One arcminute corresponds to about 2.3xl0®m at 
the comet. One pixel corresponds to 2.45". a) Mosaic of the primary images. The 
surface brightness is in MJy/sr, and moving objects close to the comet orbit have 
been removed, b) The same image after subtraction of the shadow taken six days 
after the primary. The trail was still in the FOV when the shadow was made and is 
visible in white close to the lower margin of the mosaic. For c) and d) see caption 
of Fig. 2. For this image we used a boxcar-filter of size 5x5 pixels. 
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Fig. 4. Simulations of the 2004 WFI image with the best-fit parameters o = 0.1 and 
f^l \/Qv^ ~ Pa-iisls (a) and (b) show the brightness profile and the FWHM as 
functions of the distance from the nucleus, panels (c) to (f) show brightness profiles 
perpendicular to the projected orbit at different distances from the nucleus as indi- 
cated in the upper right corner of the plots. Separate simulations were run with only 
particles emitted at least one orbital period before the observation (long-dashed) 
and emitted during the most recent apparition of the comet (short-dashed). The 
dotted lines correspond to the sum of both simulations, and the solid lines show the 
result of including the flatfield and spatial filtering in the simulation. Crosses mark 
the observation. 
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Fig. 5. Simulations of the 2005 Spitzer image with the best-fit parameters a = 0.1 
and /v/ VQpr = 1-6. For details see Fig. 4. The solid lines show the result including 
the PSF and spatial filtering in the simulation. 
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Fig. 6. Simulations of the 2006 Spitzer image with the best-fit parameters a = 0.1 
and /v/ VQpr = 1-6. For details see Fig. 4. The solid lines show the result including 
the PSF and spatial filtering in the simulation. 
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Fig. 7. Geometric albedo of the dust at visible wavelengths as a function of the 
param eter Z = (e/(l — A)) in the range of 0.4 < Z < 0.8 (ISykes and Walker 
}l993)). The function p{Z) is derived from the condition that the visible and infrared 
observations must be reproduced with the same dust production rate. We plot p{Z) 
separately for the 2005 and 2006 observation, and for two extreme values of the 
phase integral q = A/p. 
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Fig. 8. Minimum and maximmn radiation pressure coefficient of particles in the 
FOVs of our observations as functions of the emission time for /v/ ^/Qpr = 1-6. 
The upper panel shows the maximum /3 of particles remaining in the FOVs of 
our observations for emission between 1959 and 1997. The lower panel shows the 
minimum and maximum /3 of particles from the apparition immediately preceding 
our observations (2002/03). The 2005 maximum curve (not plotted in the upper 
panel) is between the 2004 and 2006 curves, and the minimum curve is the same 
during all apparitions. 
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Fig. 9. Minimum Afp parameter and mass production rates for the best-fitting sim- 
ulation parameters (a = 0.1 and fv/\/Qpr = 1-6) and particles in tlie /3 range of 
2004 as shown in Fig. 8. The Afp was measured at two nuclear distances (p = 
10000 km and p = 20000 km) in simulated images at various heliocentric distances. 
The plotted values are assuming a constant phase function. For compa rison we plot 



observed values of A fp: cros s es are data from the 1982/83 apparition (jStorrs et al 



(|l992l ): lFeldman et all (|20o3 l: ISchleicherl ^200(i )). open circles are data from 1996/97 



data 



Weiler et al. (2004): 


Schleicher 


(2006)). triangles 


(20031): Schulz et al. ( 


2004a): 
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plotted for the two extreme values of the albedo shown in Fig. 7. 
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Fig. 10. Ranges of the radiation pressure coefficient (3 of trail particles as a function 
of the distance from the nucleus. The /3 ranges were evaluated at multiples of 0.25° 
in mean anomaly behind the nucleus in simulations of all three observations with the 
best-fitting set of parameters (a = 0.1 and /v/-\/ Qpr = 1-6). The solid lines refer to 
particles stemming from the apparition immediately preceding the observations (in 
2002/03), the dashed lines to particles older than one orbital period (emitted before 
2002). Both were evaluated at the same nucleus distance but plotted slightly off-set 
from each other for better visibility. The lines indicate the presence of particles hav- 
ing a given f3 at the given nucleus distance, not their abundance. The contribution 
of particles with /3 > 0.01 to the total cross-section in the 2004 simulation is on the 
order of 10%. 
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Fig. 11. Comparison of published models of the dust emission from comet CG de- 
rived from tail or trail images. Where two or more models use the same value, the 
lines are plotted with a small offset to improve visibility, a) Exponent of the dif- 
ferential size dist ri bution, b) Min i murn and maximum /3 as a function of time. For 



Fulle et all (120041 ^ : iMoreno et al.l (l2004l'l and this paper, the /3 r ange used to fit the 
images is plotted. Ilshigur d(l2QQ8|) and lKellev et"^ Eol) do not give details 

on t he time deperidence o f the (3 ra i ige. c ) Emission spe e ds. For non- isotropic mod- 
els dMoreno et all tooi ): llshigurol ^OM); IXellev et al.l toOfj . hm^ ) the speed at 
the subsolar poi r it is given, d) Mass production rate in the /3 range shown in panel 
bV lKellev et"lll (|2nn8l . liooil l use a production rate proportional to r^^ ^'^ but do not 
mention absolute values. 
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Fig. 12. Example for a pixel (hatched area) containing part of the edges of a shell 
(area covered by shell is grey). The Zq^ are defined by Zq^ = z{\ cos(j)\ = ^,yj)- The 
integration limits for this pixel are as follows: on the upper pixel boundary, (j) = 
(t>nm{y2, z) for z E [2^1,-^0 ]) while (p = <pan{y2) for the rest. On the lower pixel 
boundary, </> = (pnmiyi, z) for z G and (j) = (/>an(yi) for z G [zi,z^^^]. For 

z E [zq ,Z2], neither boundary contributes. Hence this section does not have to be 
considered at all. This effectively corresponds to a shortening of the pixel, indicated 
by the dashed vertical line. 
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Fig. 13. Different possibilities of pixel boundaries intersecting tlie edge of a shell. 
The area covered by the shell is grey. The upper and lower rows show situations 
for the upper and lower pixel boundaries, respectively. Different "cases" occur in- 
dependently for the upper and lower pixel boundary. 1: Both corners are covered 
by the shell. 2: Both corners are outside the shell. If this scenario is given on both 
pixel boundaries, the pixel is either empty (01 = (j)2) or the shell is thinner than 
the height of the pixel ((/)i/2 = and 02/i = 0)- This situation occurs for shells 
that are collapsed to a neckline. 3: The right corner is inside, the left corner outside 
the shell. 4: The inverse of case 3. 5: Both corners are outside the shell, but there 
are two shell intersections in between. The possible sixth case (inverse of 5) cannot 
occur because shells are convex. In all cases, if there is only one intersection, the 
intersection point is labelled Zq^^ on the lower pixel boundary and Zq'^ on the upper 
one. If there is a second intersection, the leftmost intersection point is labelled Zq ' 
on the lower and z^^^ on the upper boundary. 
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